Astronomy & Astrophysics manuscript no. 18 127 


©ESO 2012 


January 9, 2012 





Probing the turbulent mixing strength in protoplanetary disks 
across the stellar mass range: no significant variations. 

Gijs D. Mulders 12 and Carsten Dominik 1 ' 3 



1 Astronomical Institute "Anton Pannekoek", University of Amsterdam, PO Box 94249, 1090 GE Amsterdam, The Netherlands 

2 SRON Netherlands Institute for Space Research, PO Box 800, 9700 AV, Groningen, The Netherlands 

3 Department of Astrophysics/IMAPP, Radboud University Nijmegen, P.O. Box 9010 6500 GL Nijmegen The Netherlands 



(N 

o 

(N 
S3 

a 



PL, 

6 



> 
m 
in 



o 

(N 



X 



Preprint online version: January 9, 2012 



ABSTRACT 



Context. Dust settling and grain growth are the first steps in the planet-formation process in protoplanetary disks. These disks are 
observed around stars with different spectral types, and there are indications that the disks around lower mass stars are significantly 
flatter, which could indicate that they settle and evolve faster, or in a different way. 

Aims. We aim to test this assumption by modeling the median spectral energy distributions (SEDs) of three samples of protoplanetary 
disks: around Herbig stars, T Tauri stars and brown dwarfs. We focus on the turbulent mixing strength to avoid a strong observational 
bias from disk and stellar properties that depend on stellar mass. 

Methods. We generated SEDs with the radiative transfer code MCMax, using a hydrostatic disk structure and settling the dust in a 
self-consistent way with the a-prescription to probe the turbulent mixing strength. 

Results. We are able to fit all three samples with a disk with the same input parameters, scaling the inner edge to the dust evaporation 
radius and disk mass to millimeter photometry. The Herbig stars require a special treatment for the inner rim regions, while the T-Tauri 
stars require viscous heating, and the brown dwarfs lack a good estimate of the disk mass because only few millimeter detections exist. 
Conclusions. We find that the turbulent mixing strength does not vary across the stellar mass range for a fixed grain size distribution 
and gas-to-dust ratio. Regions with the same temperature have a self-similar vertical structure independent of stellar mass, but regions 
at the same distance from the central star appear more settled in disks around lower mass stars. We find a relatively low turbulent 
mixing strength of a = 10~ 4 for a standard grain size distribution, but our results are also consistent with a = 10~ 2 for a grain size 
distribution with fewer small grains or a lower gas-to-dust ratio. 

Key words, protoplanetary disks - stars: pre-main-sequence - radiative transfer - dust settling - grain growth 



1. Introduction 

Protoplanetary disks are the main sites of planet formation. 
Within them, the small sub-micron sized dust grains grow to mil- 
limeter and centimeter sizes and settle to the midplane, where 
they eventually form kilometer-sized planetesimals and proto- 
planets that are the building blocks of planetary systems. How 
the dust grows and settles depends not only on the grain size and 
gas density, but also on the turbulent mixing strength of the gas 
that mixes small grains back up into the disk atmosphere and 
makes them collide in the midplane. Without turbulent mixing, 
all disks would be flat and dust would grow slower. Additionally, 
turbulent mixing is an important driver for the viscous evolution 
of the disk, for radial mixing of dust, and it may also play an 
important role in the later stages of planet formation. 

Although turbulent mixing is fundamental for our under- 
standing of disk evolution and planet formatio n, the cause of disk 
turbu l ence is not comp letely understood (e.g. lBalbus & Hawlevl 
1998; Armitage1 l201 lh . Turbulent mixing i n disks is often im- 
pleme nted using the a prescription from Sh akura & Sunvaevl 
(Il973h and lPringld(ll98lh . but the value of the turbulent mixing 
strength a t urb is hard to predict from theory. Although it could 
be empirically derived from line-broadening with future obser- 
vator ies such as ALMA , this is challenging with current facili- 
ties dHughes et alj|201 lb . The value typically used and inferred 
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from disk models trying t o match accretion rates is a v i SC ous=0.01 
(e.g. lHartmann et aUll998h . Comparison of steady-state accre- 
tion disk models t o millimeter images yields aviscous = 0.5... 10 -4 
dlsella et alj|2009h . 

The turbulent mixing strength manifests itself in the degree 
of dust settling. When dust grains grow to larger sizes, they de- 
couple from the gas and start settling toward the midplane. At 
what size and height they decouple depends not only on the 
gas density and temperature but also on the turbulent mixing 
stren gth: a lower turbulent mixing strength leads to flatter disks 
(e.g.[Dullemond & Dominik 2004b). As dust decouples, its rel- 
ative velocity also increases, which leads to fragmentation, pro- 
viding a natural stop to grain growth and replenishing the small 
dust grain population that is th en mixed up into the disk surface 
dPullemond & Dominikll2005l) . 

There is by now a wealth of observational evidence support- 
ing grain growth and dust settling in protoplanetary disks. Most 
notably, mid-infrared spectroscopic observations have shown 
that grains in the wa rm surface layers have g rown beyond the 
size of a micron (e.g. lHenning & Meeusll2009l) . In addition, the 
spectral index at (sub)millimeter and centimeter wavelengths 
shows that gra ins grow up to centimeter sizes in the cold disk 
midp lane (e.g. lBeckwith & Sargenlll991l : lMannTn gs & Em ersonl 
1994). This picture ha s been confirmed with h ydr ostatic radia- 
tive tr ansfer models bv lD'Alessio et al] d2006t) and lFurlan et al.l 
(2005) for T Tauri stars, who require a depletion of small dust 
grains at the disk surface and a population of big grains in the 
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disk midplane. Without grain growth and dust settling, hydro- 
static disk models overpredict observed far-infrared fluxes. 

Interestingly enough, this seems not to be the case for Herbig 
Ae and Be stars, the intermediate mass counterparts of T Tauri 
stars. Hydrostatic disk models do well in explaining far-infrared 
fluxes without the need to deplete the upper layers of the disk 
(iDominik et al.ll2003b IDullemond & Dominikl2004ilAcke et ail 
2009), especially in Meeus group I source^]. However, the (sub) 
millimeter slopes show that a populat ion of millimete r-sized 
grains exist in these disks as well (e.g. lAcke et al]l200l) . indi- 
cating that the grain-size population is not that different from T 
Tauri stars. This opens up the possibility that these disks are less 
settled, and that disk properties - in particular the turbulent mix- 
ing strength - could vary strongly as a function of stellar mass. 
The recent discovery that disks around brown dwarfs appear to 
be more settled than T Tauri stars f its well into this picture (e.g. 
lLada et alj|2006t ISzucs et alj|20loh . although this is not always 
found dFurlan et al.l201 II) . and indicates that the observed degree 
of settling could be inversely proportional to stellar mass. 

However, the degree of settling does not directly reflect the 
turbulent mixing strength. Disks around T Tau ri stars, Herbig 
stars a nd brown dwarfs vary in disk mass (e .g. iBeckwith et alJ 
1990; IScholz et alJ 1200(3: IVorobvovl [2009) and stellar mass, 
and observations at the same wavelength probe different ra- 
dial regions in the disk because of the changing stellar lumi- 
nosity. Comparing the degree of settling using a fixed height 
or a lower scale height for bigger grains (e.g. ID'Alessio et alJ 
120061: ISauter & Wold 1201 ll) is therefore less suitable for com- 
parison along the stellar mass range. Radiative transfer mod- 
els that describe dust settling using turbulent mixing exist 
(iDullemond & Dominikll2004bllHasegawa & Pudritzll2010l) . but 
have not been used to compare to observations directly. 

To constrain the turbulent mixing strength across the stel- 
lar mass range, we focused primarily on median spectral en- 
ergy distributions (SEDs). We did this because measuring the de- 
gree of settling in individual disks is not straightf orward because 
of a number of degeneracies in disk modeling (ISauter & Wold 
and requires detailed multiwavelength modeling, includ- 
ing scattered-light images dPinte et al. 2008), bu t conclusion s 
can be drawn from larger samples dFurlan etal] 120051 120 06). 
They are available for both T Tauri stars (e.g. iFurlan et all 2005; 
ID'Alessio etall 120061) . and brown dwarfs dSzucs et all 12010). 
and a substantial sample is available to construct them for Herbig 
stars as well (e.g. lAcke et alj |2009; Juhasz et al. 201QI). 

We will describe how dust settling is implemented in our 
radiative transfer code in section [2] and compare the available 
approaches that use the turbulent mixing strength. In section [3] 
we will use this model to constrain the turbulent mixing strength 
from the median SEDs of T Tauri stars, Herbig stars and brown 
dwarfs, which were constructed in appendix lAl and IB"1 We will 
discuss what the median disk model looks like in section [4] and 
whether disks around lower mass stars are flatter or more settled 
than their high-mass counterparts. We summarize our findings in 
section |5] 
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Fig. 1. Artist impression of vertical structure iterations in 
MCMax to determine the vertical structure. Model parameters 
are displayed in black on white, parts of the code as blue blocks. 
T refers to the temperature, rho to the density. 




combined with the method of continuous absorption by iLucvl 
dl999l) . The code has been benchmarked against o ther radia- 
tive tr ansfer codes for modeling protoplanetary disks (Pinte et al. 
2009). The radial grid around the inner radius and disk wall is re- 
fined to sample the optical depth logarithmically. The SED and 
images are calculated by integrating the formal solution to the 
equation of radiative transfer by ray- tracing. 

In addition to radiative transfer, the code can explicitly solve 
for the vertical structure of the disk. With the implicit assump- 
tion that the gas temperature is set by the dust temperature, the 
vertical structure of the gas can be calculated by solving the 
equation of hydrostatic equilibrium. On top of that, the dust 
can be settled with respect to the gas, which we will describe 
in section l2~2l After calculating the new structure for the dust, 
the radiative transfer code can be run to obtain a new dust tem- 
perature, and the whole procedure (radiative transfer, gas struc- 
ture, dust settling) can be iterated until convergence is reached 
(Figure [T). The temperature structure usually converges within 
five iterations. The exact number of iterations depends on disks 
parameters. In general, very settled or 'flat' disks take longer to 
converge because they deviate more from the initial guess. 

2.2. Dust settling 

To settle t he dust with respect to the ga s, we used the for- 
malism by IDullemond & Dominikl d2004bl) . We will briefly de- 
scribe the basic assumptions, and refer to the previously men- 
tioned paper for details. W e compare our approach to that of 
lHasegawa & P udritz (20 Tob in section [2.2.21 who used the mid- 
plane formalism by iDubrulle et al.l d 19951) in combination with 
radiative transfer models. 



2. Model 

2.1. Disk model 

The disk model used in this paper is MCMax dMin et al.| [2009). 
a 2D Monte Carlo radiative transfer code. It is based on the im - 
mediate re-emission procedure from Biork man & Wood! d2001l) . 

1 We note that for group II sources, some settling may be required to 
settle the outer disk into the inner disk shadow. 



2.2.1. Theory 

The vertical motion of dust particles within the disk is a bal- 
ance of the downward motion caused by the gravity of the star 
versus the mainly upward motion caused by turbulent mixing. 
Comparing the timescales for both processes (the friction time 
ffric versus the eddy turn-over time f e dd) results in an expression 
for the Stokes number St, which describes how well the dust is 
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coupled to the gas at every location in the disk: 



St - — - - m ^ k 



tedd 4 CT p gas C s ' 



(1) 



where m/cr is the mass-to-surface ratio of the dust grain, = 
^GM*/r 3 the Keplerian frequency, p gas is the local gas den- 
sity and c s = kT/pm pmton the local sound speed for the mean 
molecular weight p. A Stokes number smaller than one means 
the dust is well coupled to the gas, whereas particles decou- 
ple from the gas at larger Stokes numbers. Using this expres- 
sion for the dust-gas coupling a nd the standard a prescription 
from lShakura & Sunvaevl d 19731) in the form of lPringld d 19811) 
(v = ac s H p ), the diffusion coefficient for a dust particle at every 
location in the disk can be written as^ 



D = awb 



c s H P 
1 +St 2 



(2) 



where H v = c s /Qk is the local pressure scale height and a tur b is 
the turbulent mixing strength. 

Taking the gas density and temperature structure from the 
radiative transfer and hydrostatic structure calculation, we can 
now solve for the vertical structure of a dust grain of certain size 
p a at every location in the disk by solving a vertical diffusion 
equation: 



dpa(z) d 



d_ (pM\\ _ d_ 
dz \ Pgas // dz 



(PaV S ett,a(z)) , (3) 



where v set t, a (z) is the local settling speed for that dust grain in 
the absence of turbulent mixing. This equation is then solved 
in a time-dependent way using implicit integration. Because 
the timescales on which the dis k settles toward equilibrium ar e 
shorter than one million years (Dulle mond & Do minik 2004b), 
we jump toward equilibrium in ten relatively large time steps. 

We assumed that the turbulent mixing strength is constant 
throughout the disk. There is no a priori reason to do so, nor 
a computational one, and there are indications for variations in 
both the ve rtical (|Fromang"&^ Nelson 2009; Sim on et al .1 1201 11) 
and radial dlsella et alJ l2009h direction. However, because we 
cannot constrain these variations from the data used in this paper, 
we kept a constant. 




Fig. 2. Dust-to-gas ratio versus relative height (z/r) in the mid- 
infrared emitting region of the T Tauri star at 10 AU using the 
self-consistent settling (top) and the midplane settling approach 
(bottom) with the same gas density distribution. The black solid 
line represents the total dust-to-gas ratio, the gray lines denote 
the individual contributions of grains of different sizes, from 
0.02 micron (dotted) to 6 mm (solid). The dotted line represents 
a fully mixed disk, with a dust-to-gas ratio of 0.01. The radial 
t=1 surface in the optical is located at a relative height of ap- 
proximately 0. 16 and 0.2. 



2.2.2. Impact on disk structure 

The main dif f erence between this approach and that by 
Dubr ulle et alJ d 19951) is that we used the local sound speed 
and gas density to calculate the dust-gas coupling. Their ap- 
proach is an analytical solution for modeling the dust distri- 
bution in the isothermal midplane, and uses only the midplane 
density and temperature. Hence, we will refer to the latter as 
midplan e settling, though it is som etimes applied outside this 
regime dHasegawa & Pudritzl 12010). To compare the different 
approaches, we have displayed the vertical structure at the first 
iteration in a vertical slab centered at 10 AU in Figure|2] Because 
this model has not been iterated, the vertical structure of the gas 
is calculated from the optically thin dust temperature, and corre- 
sponds to a Gaussian. 

2 W e used a Schmidt num ber of 1+St 2 followinglYoudin & Lith wick 
( 120071) . rather than 1+St from Dullem ond & Dominik! d2004bl) . Because 
this does not change the point where particles decouple at St=l, it has a 
very limited influence on our model. 



The well-mixed model has a uniform dust-to-gas ratio of 
0.01, whereas the settled models show a clear trend from an in- 
creased dust-to-gas ratio near the midplane to a decreased ratio 
near and above the surface, which is located around z/r~ 0.18 
(radial t = 1 surface in the optical). This is caused by stratifica- 
tion, the settling of different dust species to different heights in 
the disk. In the midplane, the midplane settling approach agrees 
very well with the self-consistent settling, whic h is also the re- 
gion that the former has been constructed for dDubrulle et al.l 
Il995h . 

Closer to the disk surface, the two settled models start to 
deviate. With increasing height above the midplane, the gas den- 
sity decreases. For the self-consistent settling, this means that 
smaller particles will decouple from the gas. For each particle 
size, there is a characteristic height above the midplane where 
the coupling is so weak that they are almost fully depleted 
dDullemond & DominikH2004bh . For particles down to a size of 
roughly 1 micron, this decoupling takes place below the disk sur- 
face, but smaller grains can still be found above the disk surface. 
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Fig. 3. Dust temperatures versus relative height in a vertical slab 
centered at 10 AU for a T Tauri star, showing that stratification 
in dust size leads to a stratification in dust temperature. Shown 
are the self-consistent model (solid line), the settling recipe from 
Dubrulle (dotted line) and a well-mixed model without settling 
for comparison (dashed line). 



The midplane settling approach shows a different behavior, 
because the dust-gas coupling is calculated in the midplane only. 
Particles do not decouple from the gas at a specific height, but 
keep following a Gaussian distribution all the way through the 
disk surface. This gives an almost constant dust-to-gas ratio for 
the smallest particles (< 1 jum), and a much weaker depletion for 
the intermediate sizes. 

Settling the dust also impacts the temperature structure of 
the disk (Fig.O, and therefore its density structure. Non-settled 
disk models have a more or less constant temperature distribu- 
tion below the disk surface (Fig. [3] dash-dotted line). Big grains, 
however, are efficient coolers, and settling them leads to a colder 
midplane and warmer surface. The stratification in dust sizes 
also leads to a gradual increase in temperature from the mid - 
plane to the disk surface (see also Hasega wa & Pudritz] |2010). 
Because the temperature in and above the midplane is no longer 
isothermal, it becomes important to solve for the vertical struc- 
ture of the disk using the local scale height, which increases with 
height. Using only the midplane temperature to calculate the ver- 
tical structure of the disk underestimates the scale height within 
the midplane by up to a factor of two (Fig. [3}. This leads to flat- 
ter disks with - in turn - even lower midplane temperatures (Fig. 
[3] dotted line). 

2.3. Dust grain properties 

Dust grains in our model are characterized by three different pa- 
rameters: composition, size and shape. Different grain sizes set- 
tle to different heights in the disk (See Fig.[2]i, and their opacities 
are displayed in figure |4] 

2.3.1. Composition and shape 

Our dust grains consist mainly of amorphous silicates, mea- 
sured for t he compos i tion of interstellar dust toward the Galactic 
center by [Min et all (120071) . co-added with with a continuum 
opacity source, for which we used amorphous carbon. The ra- 
tio of amorphous carbon / amorphous silicates influences the 
peak strength of the 10-micron silicate feature. We used a mass 
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Fig. 4. Dust opacities for the different grain sizes used in this 
paper. 



fraction of 15 % carbon, a value found from modelin g this ra- 
tio in Herbig Ae/Be stars dMeiier et al.1 12005L 120091) . The fi- 
nal composition, wi th reference to the op tical constants, is 
11.73 % MgFeSiOd (lDorschneretal.i l 1995b. 3 2.6 % Mg 2 Si0 4 
dHenning &Stognienkol 19961) 36.5% MgSiO ? (iDorschner et al.1 
119951). 1.5 % NaAl SioOfi dMutschke etalJ Il998h . 15% C 
(Prei bisch et al.l [T993h . The shape of our particles is irregular, 
and appro ximated using a distribution of hollow spheres (DHS, 
iMin et ai1l2005l) . We used a vacuum fraction of 0.7, though this 
choice does not influence our results because we do not focus on 
detailed feature shapes. 



2.3.2. Size 

Observations of protoplanetary disks show evidence for 
(sub)micron and millimeter-sized or larger grains. Intermediate 
sized grains are n ot observed, but are predicted by dust coagula - 
tion models (e.g. IWeidenschillingll 1 984t iBirnstiel et"aTH 20 10a). 
We therefore used a grain size distribution ranging from sub- 
micron to millimeter sizes that follows a power law, f(a) = aT q , 
and used a default index of q — 3.5 si milar to the Mathis , Rumpl 
and Nordsieck (MRN) distribution (Mathi s et aDll977h . To al- 
low settling of different grain size to a different height in the 
disk, we divided our size distribution into two bins per order of 
magnitude. We define the grain radius for the settling calculation 
with the logarithmic mean of the minimum and maximum grain 
sizes to (i.e. a = 1.78 micron for grain size bin between 1.00 and 
3.16 micron). The opacities of these grains were calculated us- 
ing 10 sizes per bin to avoid strong resonant features that could 
result from using a single grain size. 



2.4. Disk parameters 

Because our disks are in vertical hydrostatic equilibrium, the 
only free fitting parameter to control the vertical structure for a 
given mass is the turbulent mixing strength a. The radial dis- 
tribution of dust and gas is parametrized by a power-law as 
S(r) = r~ p with the commonly used p — 1 for a steady-state 
accretion disk. The inner radius is placed at a dust evaporation 
temperature of 1500 K, whereas the outer radiu s is fixed to 300 
AU in accordance with iD'Alessio et alj d2006). The total disk 
mass is a fixed fraction of 1 percent of the stellar mass (e.g. 
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Parameter 


Value 


Range 


P 


1.0 


0...2.0 


M dlsk 


0.01 M sar 


0.001. ..0.1 


dust-to-gas 


0.01 


0.01... 1 


T 

" evap 


1500 


1000.. .2500 


Rou, [AU] 


300 


50... 1000 




0.01 


0.01... 10 


amaxipm] 


1000 


10 ... 10 s 


q 


-3.5 


-2.5 ... -4.5 




0.01 


10~ 6 ... 1.0 



Table 1. Parameters of our disk model, and ranges explored. 



Parameter 


Brown Dwarf 


T Tauri 


Herbig Ae 


Teff [K] 


3000 


4000 


8500 


L, [L ] 


0.025 


0.9 


21 


M. [M ] 


0.08 


0.5 


2.0 


Rta [AU] 


0.011 


0.07 


0.30 


M dust [Ms] 


8 x 10~ 6 


5 x 10" 5 


2 x 10~ 4 


distance [pc] 


165 


140 


140 


M^ [M /yr] 





3 x 10~ 8 





T halo 








0.14 



Table 2. Stellar-dependent model parameters. 
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Fig. 5. Observed median SED of T Tauri stars in Taurus (dia- 
monds) with quartiles (gray area). Overplotted is our best-fit disk 
model with a turbulent mixing strength of a = 10" 4 and an accre- 
tion rate of 3 X 10~ 8 M /yr (solid line). Also plotted are a model 
with a = 10~ 2 (dotted line) and without accretion (dashed line). 
The gray line denotes the stellar photosphere. The inset shows 
the millimeter regime. 



IScholz et al.l2006tlWiiliams & Ciezall201 lh . assuming a dust-to- 
gas ratio of 1:100. 

3. Fits to median SEDs 

To compare our disk model to the three different data sets, we 
used median SEDs. Instead of fitting all individual sources or 
comparing them to a model grid, we only fitted our model to 
one median S EP to extract 'typical' fit parameters. Furla rTetaD 
( 20051 12006b have shown that this approach works well, taking 
into account the observed spread in mid-infrared colors. 

When fitting the data, we tried to keep as many parameters 
fixed as possible, see TableQ] Parameters that could not be fixed 
because of the different stellar properties are given in Table [2] 
These are the inner radius, which we kept at 1500 Kelvin and 
the dust mass, which is a fixed fraction of the stellar mass of 
10~ 4 rvL. 

3.1. T Tauri stars 

The Taurus m edian SED has been fitted using a hydrostatic 
disk model bv lD'Alessio et all {2006), which demonstrated the 
need for grain growth and dust settling. Our goal is to re- 
produce these results, but now in the context of a more self- 
consistent settling approach to constrain the turbulent mixing 
strength, and not an arbitrary height for the midplane layer 
of big grains. We have extended the existing data set from 
Furlanetal. (2006) with sub-millimeter fluxes at 850 yum us- 
ing the catalog f r om | Andrews & Williamsl ([2005 ). The data from 
ID' Alessio et all (Q999) include some non-detections and upper 
limits (see appendix lAt. The updated millimeter median is well 
fitted by a dust mass of 5 x 10~ 5 M , consistent with the result 
from lAndrews & Williamsl (120051) . 

To fit the median SED in the near-infrared, we needed to in- 
clude viscous heating (Pig.[5j. We did not self-consistently solve 
the radial structure of the disk, but calculated the viscous heat- 
ing with a mass ac cretion rate of 3 x 10~ 8 M /yr, the same as 
Furl anetal.1 12006). To fit the mid-infrared flux, we flattened the 



disk by settling the dust. A turbulent mixing strength of a — 10~ 2 
- a value commonly used for viscous accretion disks - over- 
predicts the mid-infrared fluxes longward of 20 micron, which 
arises from the region around ~10 AU. We need to lower the 
turbulent mixing strength to achieve a good fit, which results in 
a = 1Q- 4 . 

Because the 20 micron flux mainly probes the height of the 
disk surface, there are also other ways to lower it, see also sec- 
tion [43] We can reduce the number of small grains, which will 
lower the optical depth in the upper layers and decrease the visi- 
ble disk surface. Decreasing the total disk mass would be incon- 
sistent with the millimeter flux, but we can modify the grain size 
distribution because it is dominated by big grains. Decreasing 
the power-law index of the grain size distribution removes small 
grains from the disk surface, but does not significantly affect the 
mass in big grains near the disk midplane. Changing the power- 
law index to 3.25 is sufficient to fit the median SED with a tur- 
bulent mixing strength of a = 10 -2 . 

3.2. Herbig stars 

Spectral energy distributions for a large number of Herbig stars 
are available, but a median SED has not been constructed before. 
We constructed one in appendix[Bl which we present in figure[6] 
To test the assumption that more massive stars have less settled 
disks, we retained as many disks parameters from the T Tauri 
stars as possible. We changed (see also table [2]l the disk mass - 
which scales with the stellar mass and is increased with a factor 5 
with respect to the T Tauri star, consistent with the submillimeter 
data - and the inner radius - which also increases to be consistent 
with dust evaporation. 

Although accretion does not play a role in the SEDs 
of Herbig stars because of the increased stellar luminos- 
ity and larger inner radius, the near-infrared fluxes (2-8 mi- 
cron) of hydrostatic d isk models underpredict the observations 
(iVinkovic et al.l I2006L see also Fig. [6j. Different mechanisms 
have been proposed to explain the difference, from an in- 
creased rim scale height (IVerhoeffetalJl2010t lAcke et al.l l2009) 
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Fig. 6. Observed median SED of Herbig stars (diamonds) with 
quartiles (gray area). Overplotted is our best-fit disk model with 
a turbulent mixing strength of a — 1CT 4 and a dust shell with 
optical depth t = 0.14 (solid line). Also plotted are a model with 
a = 10~ 2 (dotted line) and without a dust shell (dashed line). 
The gray line denotes the stellar photosphere. The inset shows 
the millimeter regime. 
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Fig. 7. Observed median SED of brown dwarfs in Chamaeleon 
I (diamonds), where the 24 micron and millimeter point should 
be treated as upper limits. Overplotted is our best-fit disk model 
with a turbulent mixing strength of a = 1CT 4 (solid line). Also 
plotted for comparison is a model with a = 10~ 2 (dotted line). 
The dashed line denotes the stellar photosphere. The inset shows 
the millimeter regime. 



to halos (IVinkovic et al l l2006t IVerhoeff etaT1l201ll) to mate- 
rial within the dust evapo ration radius (e.g. iKraus et alj [2008 ; 
iTannirkulam et al.1 l2008al) . Although an artificiall y increased 
scale height works well for the mid-infrared SED (Ack e et aD 
2009), it does not work well when fitting simultaneously the 
near-infrared region of the Herbig median SED (see appendix 
P . Because an inner gaseous disk is beyond the current capa- 
bilities of our code, we modeled the near-infrared flux deficit of 
the hydrostatic model with an optically thin spherical halo. We 
emphasize that this halo is a parametrization of optically thin 
emission from the inner regions, including dust or gas within the 
inner rim, and refer the interested reader to the appendix. 

The optically thin halo plays the role of accretion in T Tauri 
stars: it provides near-infrared flux without affecting the flux 
longward of 20 micron, in contrast to the puffed-up inner rim. 
We used a small halo of 0.3... 1 micron grains between 0.14 and 
0.3 AL0 that has an optical depth of 0. 14 in the visual and a dust 
mass of 4 x 10~ 12 M (Figure|6]l. As a result of the fitting proce- 
dure, we found the same turbulent mixing strength for the disk 
of a — 10~ 4 as for the T Tauri stars. A higher mixing strength 
of a — 10~ 2 overpredicts the Spitzer IRS fluxes from 20 to 40 
micron and the IRAS flux at 60 micron. 



3.3. Brown dwarfs 

Disks around brown dwarfs are more difficult to detect owing 
to their low luminosity. Complete samples down to photospheric 
value s only exist up to 8 micron, and are incomplete at 24 micro n 
(e.g. lLada et al.ll2006l IScholz et al.ll2007l iLuhmanet al.l 12 008). 
Furthermore, very few m illimeter detections exist dKlein et al.l 
l2003HSchorz et all2 006). We used the median SED constructed 
by ISzucs et alj d2010l) for Chameleon I, where the median 24 
micron flux should be treated as an upper lim it. For convenience , 
we also plot the millimeter upper limit from Klei n et al . (2003), 
scaled to the distance of the Chameleon I star-forming region. 



3 Because the halo spans a very narrow temperature range, its near- 
infrared spectrum is virtually indistinguishable from that of a puffed-up 
inner rim (see Fig. lC. It . 



Unlike the Herbig and T Tauri stars, hydrostatic disk mod- 
els of brown dwarfs do not show a flux deficit in the near- 
infrared compared to the observed median. We scaled the disk 
mass and inner radius to the brown dwarf regime (see Table 
The ste llar photosphere was f itted by a NextGen stellar atmo- 
sphere (Hauschildt et al. 1999), and we adjusted the luminosity 
down to 0.025 L Q . We fitted the median shortward of 8 micron 
without any additional adjustments (see figure|7]). The assumed 
dust mass of 8 x 10~ 6 M Q is consistent with the upper limit 
from Klei n et al.l d2003l) . and amounts to only 2-3 earth masses 
of solid material. 

We found that a turbulent mixing strength of a = 10~ 4 is con- 
sistent with the 24 micron upper limit, whereas a mixing strength 
of a — 10~~ 2 clearly overpredicts this upper limit. We found the 
same turbulent mixing strength as for the earlier spectral types. 



4. Discussion 

Although we found the same turbulent mixing strength for the 
three different samples, this does not necessarily mean that these 
disks are equally flat. The vertical structure is a product of more 
than just the turbulent mixing strength, most notably the local 
density, temperature and Keplerian frequency. In the following 
section we will describe the disk properties, and show that disks 
around lower mass stars appear mo re settled at a specific radius, 
consistent with ISzucs et all (12010). but that the structure in the 
mid-infrared emitting region is remarkably self-similar. 

We will also discuss why our choice of grain size distribution 
does not affect our conclusion that the turbulent mixing strength 
does not vary along the stellar mass range in section 1431 and its 
implication for the planet-forming potential of the disk in section 
14.41 The dependence of the fitted turbulent mixing strength on 
other disk parameters is discussed in section^ 



4.1. Disk properties at fixed radius 

To compare the dust and gas distribution in our best-fit mod- 
els, we display their pressure scale height and the radial r = 1 
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Fig. 8. Location of the dust and gas in the disk of the best-fit model of a T Tauri star (center), Herbig star (left) and a brown dwarf 
(right). Displayed in gray is the gas located within one (dark) and two (light) pressure scale heights. The solid lines marks the 
location of the radial t = 1 surface in the visual. 



surfaces in figure [8] The pressure scale height is highest in the 
brown dwarf disk and lowest in the Herbig star, reflecting the 
differences in stellar mass. Although temperatures are also lower 
for less massive stars, this is not enough to compensate the lower 
gravitational potential in which the disk resides, and the gas re- 
mains more flaring. 

The dust, on the other hand, follows a different path: the 
t = 1 surface in the outer disk reaches a relative height (z/r) of 
0.17, 0.15 and 0.14 for the Herbig staiQ, T Tauri star and brown 
dwarf, respectively. This indicates that disks around less massive 
stars are a little flatter, though the difference remains within 20%. 
Compared to the much higher pressure scale height at the same 
radius, the disk of lower mass stars therefore appears more set- 
tled, because the disk surface reaches down to one pressure scale 
height for the Herbig star, but to less than half a pressure scale 
height in the brown dwarf. This can be explained by lower (sur- 
face) densities and temperatures (see figure [TOh, [TTh and Table 
[3]) that make the dust decouple closer to the midplane, as well as 
a lower optical depth. 

The shape of the disk surface can also be viewed in 
scattered-light images. We display the radial profiles of synthetic 
scattered-light images at 0.5 yum in figurefTSk. The radial profiles 
are scaled to the photospheric flux at the same wavelength, such 
that it traces only the shape of the disk surface, and not the prop- 
erties of the central star. A clear trend is visible at all radii, with 
higher fluxes and shallower slopes for Herbig stars and lower 
fluxes and steeper slopes for brown dwarfs (see also Table [3]). 
This indicates that the disk surface at a specific location is less 
flaring. 

4.2. Self-similar solutions in the MIR emitting region 

Because the mid-infrared emission probes a different region in 
all three samples, it is also interesting to compare regions with 
similar temperatures, which are located at different radii for dif- 
ferent stars. These regions can be identified by using that the dis- 
tance r, at which a specific flux from the star is received, scales 
with the square root of the stellar luminosity L». By scaling the 
radius to the value of the stellar flux, rj = rj V£*/L , we can 
construct dimensionless equivalents of our disk properties (de- 
noted by a "), and check if they are self-similar across the stellar 
mass range. 

For example, the optically thin dust temperature for the T 
Tauri star can be described as (see Fig. |9^) 



WO = 307 K * (r/AU) 



-0.5 



4 It has to be noted that the halo increases the calculated r - 
for the Herbig star. 



(4) 



1 surface 



By scaling the radius to the intensity of the radiation field, we 
can retrieve a self-similar solution for the optically thin dust tem- 
perature profile 



TdustW = 307 K * (r T /AU * ^JLJL B ) 



-0.5 



= 315K*(r T /AU)-°- 5 = f dust (r T ), 



(5) 



where the temperature at 1 AU for the T Tauri star is close to the 
selfconsistent solution because its luminosity is close to 1 L Q . As 
can be seen in figure |9j), these scaled solutions agree to within 
20% for the T Tauri star, Herbig star and brown dwarf. 

4.2.1. Self-similar temperature and surface density 

With the same procedure, we can also compare different quan- 
tities to see if they are also self-similar in all three samples. We 
find that also the midplane temperatures for all three samples 
(see Fig.fTOb) can be described with a self-similar function: 

n-0.44 



7mid(r T ) = 90 K * (rr/AU)" 



(6) 



which is a strong indication that this must be a consequence 
of a self-similar vertical structure across the stellar mass range, 
which will we explore below in more detail. Interestingly, this 
scaling law also works for the surface density (Fig. [PTT i. such 
that regions with similar temperatures also have similar surface 
densities: 

S(r T ) = 0.25 g/cm 2 * (rr/AU)" 1 , (7) 

This is because across the three samples, the luminosity scales 
roughly with stellar mass squarecfl Since the disk mass is a fixed 
fraction of the stellar mass, the surface density scales as 



'1AU 



Mo 



-1AU 




(8) 



However, the region corresponding to the same temperature lies 
farther away from the star, compensating for the change in disk 
mass: 

-t 



S(r T ) = Huu r T ~ 



(9) 



Thus the self-similar solutions for the surface density agree to 
within 20 %. 



5 Although the main-sequence mass luminosity relation for solar 
mass stars follows L oc M 3 , the T Tauri stars in our sample are 
pre-main sequence and do not follow this relation, whereas the brown 
dwarfs are in a different scaling regime altogether. 
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Fig. 9. Radial temperature profile of the optically thin dust for the best-fit model of a T Tauri star (solid line), Herbig stars (dotted 
line) and a brown dwarf (dash-dotted line). The left panel shows the real temperature, the right panel shows the self-similar solutions. 




Fig. 10. Radial temperature profile in the midplane for the best-fit model of a T Tauri star (solid line), Herbig stars (dotted line) and 
a brown dwarf (dash-dotted line). The left panel shows the real temperature, the right panel shows the self-similar solutions. 




Fig. 11. Radial surface density profile for the best-fit model of a T Tauri star (solid line), Herbig stars (dotted line) and a brown 
dwarf (dash-dotted line). The left panel shows the real surface density, the right panel shows the self-similar solutions. 




10 100 10 100 

distance [All] R[AU]/Vl 

Fig. 12. Radial intensity profile in scattered light at 0.5 micron for the best-fit model of a T Tauri star (solid line), Herbig stars 
(dotted line) and a brown dwarf (dashed line). The left panel shows the real intensity as fraction of the stellar flux at that wavelength 
F SC at(r)/F» versus r), the right panel shows the self-similar solutions (F scat (r)/F» • (L»/L ) versus r/ VL»/L Q 
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4.2.2. Self-similar scale height 

If we now look at a vertical slab in the disk located at rx = 
rj VlT/Zq) we find that two of the main parameters for the set- 
tling calculation in equations ([TJ and (f2]i - the gas surface density 
(O and temperature (|6]l - are self-similar. The third important pa- 
rameter, the Keplerian frequency, is not the same, but scales with 
VWL 



Q k (r T ) = 




(10) 



GM G 



yjUILo = Shir) jL*/L e 



which is equivalent to t/tj. The gas pressure scale height at this 
location is therefore not self-similar, but the relative scale height 
Hp/ris: 



H p (r T ) c s (r T ) 



c.(r) 



c s (r) ff p (r) 



rr r T £2 k (r T ) r T (Q k (r) r/r T ) rQ k (r) r 

(11) 

This means that if we express the height above the midplane (z) 
in terms of the relative height (z/r), the local scale height, and 
therefore the vertical structure of the gas, is self-similar: 



H p (r T ,ZT/r T ) = H v (r T ,zj/r T ). 



(12) 



4.2.3. Self-similar dust settling 



The vertical structure of the dust in the mid-infrared emitting 
region at rj - 10AU is shown in figure [13] Like that of the gas, 
the vertical structure of the dust also appears self-similar when 
expressed in terms of its relative height z/r. This is because also 
the Stokes number (Eq. ([TJ) 



« 3 m 

St(r T ,z T /r T ) = - 



Q k (r T ) 



4 crp gas (r r ,z T />T) c s (r T ,ZT/r T ) 
3 m Q k (r)/ ^f^7^ 



(13) 



4 °~ Pgas(r, z/r)/ *4UTQ c s (r, z/r) 
= St(r,z/r) 

is self-similar because the gas density scales with l/H p . 
same is true for the diffusion coefficient (0 



The 



£>(r T ,ZT/r T ) = amrb 



a'turb" 



c s (r T , ZT/r T )tfp(r T , Z T /r T ) 

1 + St (r T ,z T /r T ) 
c s (r T , z T /r T )//p(r T , z T /r T ) 



(14) 



1 +St 2 (r T ,z T /r T ) 



= D(r,z/r) 



because it is only a product of self-similar quantities when ex- 
pressed in units of r? and Zr/rr- Consequently, the vertical struc- 
ture of both gas and dust in the mid-infrared emitting regions of 
protoplanetary disks are self-similar. 

This result is consistent with that o flFurlan et alj J20TTJ) . who 
also founnd no variation in the degree of settling betw een brown 
dwarfs and T Tauri stars, but partly contradicts ISzucs et aD 
(l20Toh . These authors also found the same relative scale height 
(H/r) for the dust in the mid-infrared emitting region. However, 
the required reduction in dust scale height compared to their 
fully flared models implies a difference in gas scale height be- 
tween brown dwarfs and T Tauri stars, which we do not find. 




0.20 



Fig. 13. Dust-to-gas ratio at rj = 10 AU for the best-fit model 
of a T Tauri star (solid line), Herbig star (dotted line) and brown 
dwarf (dashed line). The gray line denotes a well-mixed model 
with a dust-to-gas ratio of 0.01. The surface temperature at 
this location corresponds to roughly 100 K, and therefore emits 
mainly in the mid-infrared. 



Disk properties 



Parameter 


Brown 


T Tauri 


Herbig 


Self- 




dwarf 




Ae/Be 


similar 


Tiau [K] 


36 


86 


185 


90 


q 


0.45 


0.44 + 


0.43 


0.44 


Siau [gr/cm 2 ] 


0.04 


0.24 


0.95 


0.25 


Iiooau [mJy/AU 2 ] 


1.110 9 


1.3-10- 6 


4.2- 10" 4 


2.9- 10~ 5 





3.1 


2.9 


2.6 


2.9 



Table 3. Disk properties for our best-fit models. The midplane 
temperature is fitted by T mp (r)=TiAu 1 ' q - The intensity in scat- 
tered light by I(r)= Iiooau r ° outwards of 70 AU. ' The midplane 
temperature in T Tauri stars shows a clear break at 2 AU caused 
by viscous heating (fig [T0h). and at smaller radii it is best de- 
scribed as T(r)= 125 r J 093 . 



The self-similarity of vertical structure is further illustrated 
by synthetic scattered-light images that trace the shape of the 
dust disk surface (Fig. [T2l . If we plot the radial intensity pro- 
file as a function of the scaled radius rj, corrected for the rela- 
tive height (F sca t/F,, * L„/L ), we see that the radial profiles for 
the three different stars line up nicely (Figure [T2b). and that the 
shape of the disk surface is indeed self-similar. 

Since the amount of radiation reprocessed at every radius 
is set by the covering fraction (z/r) - which is the same for all 
stars - even the SED in the mid-infrared becomes self-similar 
(figure [T4l . The scaled model SEDs line up to within 20% in 
the mid-infrared at 20-30 micron where we fitted the SEDs. The 
similarity extends into the far-infrared (within 40%) and even the 
near-infrared up to 3 micron (within a factor of 2). The latter is 
surprising and probably coincidental, given the different nature 
of the emission mechanism at that wavelength: viscous heating 
for the T Tauri star, optically thin emission for the Herbig star 
and mostly photospheric flux for the brown dwarf. 



6 The relative height z/r scales with rj/r 
intensity in scattered light scales with (r-i/r) 1 ■ 



1 / y/LJL a , while the 
l/(L,/L e ). 
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Fig. 14. Best-fit model SEDs for the Taurus median (solid line), 
Herbig median (dotted line) and brown dwarf median (dashed 
line), scaled to a luminosity of 1 L Q and a distance of 140 pc. 



4.3. Grain size distributions 

Throughout this article, we have assumed that there are no radial 
gradients in the grain size distribution. Models o f grain growth 
predic t that grain sizes decrease with ra dius (e.g. | Birnstiel et al.l 
2010b), w hich has also be en observed by Guilloteau et al. (201 1) 
and|Banzatti et al. (2011). Additionally, disks around later type 
stars are colder and less massive, leading to smaller grain sizes. 
Modeling the grain size distribution is beyond th e scope of this 
paper , though a recipe for doing so is available (Birnstiel et al. 
1201 ll) . For now, we will only discuss if we expect to see differ- 
ences in the grain size distribution in the mid-infrared emitting 
region between the three different samples. 

Because the vertical structure in the mid-infrared region is 
self-similar, we investigated if this is also true for the grain- 
size distribution and its planet-forming potential. Grains in pro- 
topla netary disks are in a c oagulation-fragmentation equilib- 
rium (Weidenschilling 1984), driven by turbulence. According 
to Birnsti el et al.l d201 ll) . the most important parameters that de- 
fine the size distribution are the turbulent mixing strength, the 
gas surface density, the temperature, and three parameters that 
refer to the microphysics of dust: the fragmentation velocity Uf , 
the solid density of dust p s and the size distribution of fragments 
We assumed the last three to be the same across the three sam- 
ples. 

We find no variations in the turbulent mixing strength (sec- 
tion[3]l. Regions with similar temperature also have a similar sur- 
face density (section l4~2l . The grain size distribution should then 
also be self-similar. We can see why this is the case by looking at 
the expression for the maximum grain size for turbulence-d riven 
growtfQ - as given by equation (12) of Birn stiel et al.l d2009l) : 



K(T) = 



2(r T )Wf 



T.(r)u 2 f 



?ratui-bP s Cs(r T ) na tulb p s ci(r) 



= a„ 



kW, (15) 



which is also self-similar when looking at regions with similar 
temperature (ij = xj -\JLJLq). This means that also the grain 
size distribution in the mid-infrared emitting region is not ex- 
pected to vary across the stellar mass range. 



7 which is the relevant regime for particle sizes smaller than a mil- 
limeter 



4.4. Planet-forming potential 

If the same disk properties are found at a different location in 
the disk around heavier stars, does this mean that the resulting 
planetary systems are simply scaled versions of each other? 

To compare the later stages of planet formation, we looked 
at the Hill sphere and feeding zone of a forming planet, to see 
how much mass is ava ilable for its formation, in a way similar to 
iRavmond et alJ (120071) . The Hill sphere is given by 



nip 
3Ml 



(16) 



where nip is the mass of the planet. Expressed in the same units 
that result in self-similar solutions for the vertical structure (rx = 
r/ VL*/L ) it becomes 



^h(^t) = r 



nip 



3mj ^^LJ^ ynm 



(17) 



So the size of the Hill sphere is not self-similar, though its de- 
pendence on luminosity is weak. It decreases in size for less lu- 
minous stars, despite the lower gravity of the central star. The 
total mass available in the feeding zone deviates more because it 
scales with an extra factor 2nrj\ 



M feed (r T ) = 2nr T T.(r T ) * 2R H (r T ) 
= 2nl.(r)r/ yjLJL *2 



Rn(r) 



(18) 



L,/L-J /6 M feed (r). 



The size of the feeding zone therefore scales almost linearly with 
luminosity. This means that although the initial conditions for 
planet formation are the same across the steller mass range, this 
is not true for the later stag es of planet formation, in agreement 
with lRavmond et all d2007h . 



4.5. Dependence on disk parameters 

The method we present here is useful for constraining variations 
in the turbulent mixing strength, but it is more difficult to con- 
strain its absolute value. The mid-infrared flux with which we 
fitted our models traces the height of the disk surface, which 
can also be influenced by other factors than the turbulent mixing 
strength. However, because our solutions for the mid-infrared 
emitting region are self-similar, these uncertainties will affect 
our estimate of the turbulent mixing strength in the same way 
for all stars, and do not influence our conclusion that it does not 
vary across the stellar mass range. 

In this section we will explore how other parameters can af- 
fect our estimate of the turbulent mixing strenght. In particu- 
lar, we will explore how to keep the turbulent mixing strength 
fixed at a tm \, = 0.01 and fit the three median SEDs by vary- 
ing these other parameters. As already mentioned in section lXTl 
small changes in the grain size distribution also affect the height 
of the disk surface. Using a flatter slope for this distribution re- 
moves small grains from the disk surface, and we can fit the 
median SEDs with a power-law index of 3.25. Another param- 
eter that directly influences the SED in the mid infrared is the 
gas-to-dust ratio. A lower ratio leads to a weaker coupling of the 
dust and hence to flatter disks. We fitted a gas-to-dust ratio of 1 
for a tur b = 0.01. 

One parameter that we cannot change is the dust mass, be- 
cause this would be inconsistent with the millimeter photometry. 
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As a consequence, all parameters that influence the millimeter 
opacity require a refitting of the dust mass. For example, we can 
reduce the opacity in the upper layers by reducing the carbon 
content of our dust grains, but this also reduces the opacity at 
millimeter wavelengths. If we then increase the dust mass to fit 
the flux at 850 /mi, we again increase the opacity in the upper 
layers. The two effects cancel out, and there is no net decrease in 
the mid-infrared flux. In a similar fashion, the maximum grain 
size a max does not influence the estimate of a tur b because it de- 
creases the opacity in the upper layers and at millimeter wave- 
lengths simultaneously. 

The same is true for parameters that influence the distribu- 
tion of mass in the outer disk: the outer radius and surface den- 
sity powerlaw. Although a smaller disk has a smaller covering 
fraction for the same surface density, a fit to the millimeter data 
requires a higher surface density, compensating the decrease in 
covering fraction. The same is true for steeper surface density 
powerlaws. 

4.6. Spread in disk parameters 

iFurlan et al.l ((2005) have shown that the spread in observed disk 
colors in the mid infrared is well reproduced by a spread in incli- 
nation, accretion and surface layer depletion. However, millime- 
ter observatio ns show that there is also a considerable spread 
in disk mass dAndrews & Williamsl l2005). which will affect the 
amount of mid-infrared emission through changing dust surface 
densities and stronger or weaker dust-gas coupling. We investi- 
gated if this is sufficient to explain the observed spread in disk 
colors, without the need for a spread in surface layer depletion. 

We ran a series of disk models based on the best fit to the 
Taurus median SED with or tur b = 10~ 4 , varying only the disk 
mass over two orders of magnitude. The result can be found in 
figure[T5] Although the spread of this grid reproduces the region 
between the lower and upper quartile at 20-40 micron well, it is 
much too large for the quartiles in the millimeter. This means 
that the real spread in disk mass must be smaller, and therefore 
cannot be solely responsible for the spread in observed colors. 
Albeit slightly smaller than previously estimated, a spread in sur- 
face layer depletion must still be presenfl 

This result also serves to validate one of the key assumptions 
of this paper: that a model with median fit parameters is a good 
representation of the median SED. Because the model SEDs do 
not overlap (which is also the case when varying most other disk 
parameters) the median SED is exactly equal to the SED with 
the median fit parameters. 

5. Conclusions 

We have investigated the turbulent mixing strength in protoplan- 
etary disks across the stellar mass range by comparing radiative 
transfer models with dust settling to median SEDs of three sam- 
ples of T Tauri stars, Herbig stars and brown dwarfs. Our con- 
clusions are: 

- We find no significant variations of the turbulent mixing 
strength across the stellar mass range from Herbig stars down 
to brown dwarfs. 

- We find a relatively low turbulent mixing strength of orturb = 
10 4 for an MRN-like grain size distribution up to a millime- 
ter and a gas-to-dust ratio of 100, whereas a size distribution 

8 Which would in our model translate into a spread in turbulent mix- 
ing strength or grain size distribution. 
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Fig. 15. Observed median SED of T Tauri stars in Taurus (dia- 
monds) with quartiles (gray area). Overplotted is a series of disk 
models with varying disk mass, centered around the best fit to the 
median SED. The dust mass ranges from Md us t = 5 x 10~ 4 M Q 
at the top to M dust = 5 x 10~ 6 M at the bottom, with fixed 
gas-to-dust ratio. The gray line denotes the stellar photosphere. 
The inset shows the millimeter regime. 

with fewer small grains or a gas-to-dust ratio of 1 is required 
for a higher mixing strength of a tul -b = 10~ 2 . 

- The mid-infrared emitting regions of protoplanetary disks 
are remarkably self-similar, with a similar surface density, 
vertical structure of gas and dust, and even grain size distri- 
bution. This provides comparable environments for the first 
stages of planet formation, though the later phases domi- 
nated by the planet's Hill sphere are likely different. 

- When viewed at the same distance from the central star, disks 
around later type stars appear more settled. They have a flat- 
ter disk surface as traced in scattered light. However, their 
gas scale height is much higher. 

- We have constructed a median SED for Herbig stars, analo- 
gous to its Taurus and Chameleon I counterparts. This me- 
dian SED points to high inner-rim temperatures of ~1700 K, 
and is only poorly fitted with a puffed-up inner rim, but is 
more consistent with emission from optically thin material 
from above the rim or from within it. 
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Appendix A: Taurus median 

As mentioned in section [3~Tl we have added the 850 mic r on flux 
point to the existing Taurus median from Furl an et aD (12 006). 
Althou gh the original Taurus median from iD'Alessio et alj 
(1 19991) is based on the same sources and includes millime- 
ter photometry, its coverage is incomplete. It is strongly bi- 
ased toward brighter - and thus heavier - disks and was 'meant 
to be indicative rather than definitive'. lAndrews & Will iams 
(2005) have created a nearly complete census of the Taurus star- 
forming region at millimeter wavelengths, allowing us to com- 
plete the median fluxes at 850 micron. The median flux we de- 
rive is log(vF v [erg/s/cm 2 ]) = -12.70, with the upper and lower 
quartiles at -13.05 and -12.13, respectively. These fluxes are 
roughly an order of magnitude lower than those included in the 
original median, but are consistent with the av erage disk mass in 
Taurus found by|Xndrews & Williams] (120051) . 

Appendix B: Herbig median 

The H erbi g median SEP is ba sed on the samples of I Acke et al.l 
(2009) and Juhasz et al . (2010). Because knowing the dust mass 
is crucial in measuring the degree of settling in the outer 
disk, we select ed only those sou rces with a millimeter detec- 
tion. Following lAcke et al.l (12009). we also excluded transitional 
disks, since their mid-infrared SEDs should be explained in the 
context of gaps and inner holes, rather than dust settling. 

The sample consists of 32 sources, and has a complete wave- 
length coverage in bands, B, V, J, H, K, L, Spitzer IRS long- 
ward of 10 micron, IRAS 12, 25 and 60 and millimeter wave- 
lengths. U- and M-band photometry were also included because 
they lacked only one and two measurements, respectively. The 
Spitzer spectra were reduced to narrow-band ph otometry follow- 
ing the same method as in iFurlan et al.1 (120061) . At the shortest 
wavelengths (5.7, 7.1, 8.0, 9.2 yum) the coverage is sometimes 
incomplete (lacking 5,5,3 and 3 measurements resectively), and 
was supplemented with ISO data where available. We only in- 
cluded the IRAS 60 micron fluxes, because the 12 and 25 micron 
regions is covered by Spitzer, and results in the same median 
fluxes. Although all sources have millimeter detections, not all 
are measured at the same wavelength. We therefore constructed 
one photometric point at 850 //m, and interpolated fluxes accord- 
ing to vF v (mm) = vF v (obs) * (A b s ,mm) -4 , averaging if more than 
one measurement was available. 

The median star in our sample is a Herbig A6 star, with a me- 
dian effective temperature of 8500 K, luminosity of 21 L Q and 
mass of 2 M Q . Because these stars are not in a single star-forming 
region, we scaled their SEDs to a distance of 140 parsec and the 
luminosity to the median luminosity before constructing the me- 
dian. The median including quartiles are given in table IB. ll and 
are displayed in Figsl6land lC.ll This approach results in an SED 
that has a small spread in photometry at stellar wavelengths, and 
a similar spread as the Taurus median SED at disk wavelengths. 
Typical features like the 2 micron bump and silicate features are 
clearly visible. 

Appendix C: The inner regions of Herbig stars 

As mentioned in section 13.21 the inner regions of Herbig stars 
deserve special attention in disk modeling. The strong 2 micron 
bump is a promine nt feature in most Herbig Ae and Be stars 
dMeeus et al.l l2~001). and has been explained as the fully illumi- 
nated surface of a disk truncated at the dust evaporation radius 
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band wavelength median vF„ lower quartile upper quartile 
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Table B.l. Median SED for Herbig stars, normalized to a dis- 
tance of 140 parsec and a luminosity of 21 L . 
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Fig. C.l. Observed median SED of Herbig stars (diamonds) with 
quartiles (gray area). Overplotted is a disk model with a puffed- 
up inner rim with a scale height 2.5 times the pressure scale 
height and turbulent mixing strength of a = 10 -2 (solid line). 
Also plotted are a model with the same inner rim and a modi- 
fied grain size distribution with a power-law index of 4.0 (dotted 
line) and a model without a puffed up inner rim, but with mil- 
limeter sized grains within the dust evaporation radius (dashed 
line). The gray line denotes the stellar photosphere. The inset 
shows the millimeter regime. 



bv lDullemond et all (120011) . Because of the large fraction of re- 
processed light in the near-infrared, this would require the rim to 
be 'puffed up ' , a natural cause of its high temperature down to 
the midplane (Natta et al. 200 jb. 

However, Vinkovic et al ] (I2006h showed that a rim in hy- 
drostatic equilibrium does not puff up far enough to explain 
the near-infrared fluxes. To fit the SED, the inner rim has to 
be puffed up beyond hydrostatic equ ilibrium by a factor 2-3 
dVerhoeff et al.ll2010fclAcke et alj2009h . Other ex planations have 
also been proposed, ran ging from dust halos (Vin kovic et alJ 
2006; Verhoeff et al. 1201 ll) to optically thin gas or dust w ithin 
the inner rim (Kraus et aLl l2008t iTannirku lam et alj|2008af) . as 
well as theoret ical motivation s for an increased scale height in 
the inner rim dThi et alJl20"TTh or the presence of such a halo 
dKriit & Dominikll2011l) . For the context of this paper, it is im- 
portant whether or not the added emission in the inner disk in- 
fluences the geometry and emission of the outer disk, and could 
affect our estimate of the turbulent mixing strength. A puffed-up 
inner rim casts a sha dow on the outer d isk, decreasing the flux at 
longer wavelengths dAcke et alJl2009l) . whereas material within 
the in ner rim or a halo does not cast a shadow (Mulders et alJ 
I2010h . 

Figure IC.ll shows these alternative inner disk geometries. 
When we increased the scale height in the inner rim by a fac- 
tor 2.5 to fit the near infrared part of the SED (Fig IC.ll solid 
line), it cast a shadow over almost the entire outer disk, and we 
were unable to fit the far-infrared flux even with a turbulent mix- 
ing strength of a — 10~ 2 or higher. We were able fit the far in- 
frared flux by increasing the amount of small grains, which we 
did by increasing the slope of the grain size distribution to q=4.0 
(dotted line). However, the strong effect of the inner rim shadow 
remains, supressing the mid-infrared flux. We tried several pre- 
scriptions to modify the inner rim structure, such as rounding it 
off, but they did not provide enough mid infrared flux. The rea- 
son why the puffing up the inner rim does not work well in our 



model - in contrast to e.g. lAcke etail d2009l) and lThi et all (1201 ll) 
- is that the median SED requires additional flux at 2 micron - 
whereas previous work focused mainly on the 8 micron region. 
This requires a hot inner rim, 1500-1700 K, which is consistent 
with observations and dust evaporation, but which does not pro- 
duce enough MIR emission. 

Alternatively, dust or gas within the inner rim can also con- 
tribute to the near-infrared flux. Additional components origi- 
nating in that region hav e been observed for a number of Herbi g 
stars, such as MWC 147 dKraus et alJ2008tlBagnori et all2 010). 
MWC 2 75 and AB Aurigae (ITannirkulam et alJ I2008allb1), HP 
163296 dBenistv et alJl2010h and HR599 dBenistv et alll201 ll) . 
Material within the inner rim can increase the n ear-infrared 
flux i f it generates extra luminosity from accretion dKraus et al.1 
2008) or if it captures additional radiation that would otherwise 
not reach the disk. 

For the second option, the material needs to extend suf- 
ficiently close to the star, such that starlight is captured that 
crosses the midplane and would normally not be reprocessed by 
the disk. We modeled this by extending the largest grain compo- 
nent in our model down to two stellar radii, or 0.02 AU (Fig lC.ll 
dashed line). Because even the largest grains in our model at this 
location become extremely hot, we were unable to fit the shape 
of the near-infrared SED this way, but it does capture roughly 
sufficient energy at shorter wavelengths. Because this solution 
to the near-infrared flux problem is similar to a halo in the sense 
that it captures additional energy, we did not explore it in more 
detail. We limited ourselves to modeling a halo, where we note 
that this could also represent a solution with material within the 
inner disk. 
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